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Abstract 

We present the exact solution to the moment problem for the spin- 1/2 isotropic 
antiferromagnetic XY chain with explicit forms for the moments with respect to 
the Neel state, the cumulant generating function, and the Resolvent Opera- 
tor. We verify the correctness of the Horn-Weinstein Theorems, but the ana- 
lytic structure of the generating function (e~ tH ) in the complex i-plane is quite 
different from that assumed by the i-Expansion and the Connected Moments 
Expansion due to the vanishing gap. This function has a finite radius of con- 
vergence about t = 0, and for large t has a leading descending algebraic series 
E(t) — Eq ~ At~ 2 . The Resolvent has a branch cut and essential singularity near 
the ground state energy of the form G(s)/s ~ S|s + l| _3/ ' 4 exp(C|s + l| 1 / 2 ). Con- 
sequently extrapolation strategies based on these assumptions are flawed and in 
practise we find that the CMX methods are pathological and cannot be applied, 
while numerical evidence for two of the i-expansion methods indicates a clear 
asymptotic convergence behaviour with truncation order. 
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1 Motivation 



Some novel non-perturbative formalisms based on a moment approach have been pro- 
posed recently, the t- Expansion |1| and the Connected Moments Expansion (CMX)P], 
and applied in a number of areas ranging from quantum chemistry 0, ^, ||, [5|, ||, lattice 
gauge theories f7|, |, g 0, [II], 0, [Lj to condensed matter problems 0, [OJ |T|, 0, [[§ . 



These formalisms take their starting point from a finite sequence of low order con- 
nected Hamiltonian moments, or cumulants {y n : n — 1 . . . }, and appeal to ostensibly 
exact theorems for ground state properties in order to extract appropriate averages. 
In all these applications, the sequence of cumulants was finite and relatively small, 
and the effects of truncation and the rate of convergence was completely unknown. 
This issue has particular practical importance because in many of the applications the 
convergence and final estimates were rather poor considering the effort expended in 
calculating a large number of moments. As rough rule of thumb, the initial and low 
order estimates were very good and the situation looked promising but when taken to 
higher order there was not a significant improvement, and in some cases the convergent 



behaviour seemed to be arrested [17, |19| , The expectation of an accurate, perhaps 



arbitrarily and controllably accurate, method has to date seemed to be unfounded. 
The source of these problems had not been properly identified nor any remedy pro- 
posed. Here we give the first example of an exact solution to a moment problem which 
can answer some of these questions. 

Furthermore it was not known whether these techniques were capable of yielding 
exact results in any case. This is a natural expectation since any formalism, applicable 
to general problems, must have a set of simple, but nontrivial solvable cases. There was 
another model |2(J, of quadratically coupled harmonic oscillators, for which the exact 



ground state wavefunction and density matrices are known but the moment problem 
was not exactly solved in that a general cumulant was not given. We present the 
solution to the moment problem for the isotropic, one-dimensional antiferromagnetic 
spin | XY Model. Our solution also demonstrates an alternative notion of solvability, 
from those concepts and methods usually associated with the known solvable models. 

In the t-Expansion, one constructs the cumulant generating function from its 
taylor series expansion about t = 0, 



(e^)=exp Vz, J- e/O 




where all moments are calculated with respect to some trial state, with the proviso that 
the overlap with exact ground state is nonzero. Then the Horn-Weinstein Theorem[]I[ 
states that the ground state energy and first excited state gap (within the same sector 
as the ground state) are given by the following limits 

E = lim E(t) 

t^oo 

<PE{t)/dt* (2) 
El -E =lna- dm/dt , 

where E(t) = F'(—t). Other ground state averages can be found as well by computing 
the appropriate generalised connected moments]]]] and performing the same extrapo- 
lation. This theorem was "proved" in Ref. (JT|) by assuming a simple discrete point 
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spectrum for the system, even in the thermodynamic limit, and these kinds of assump- 
tions continue to be commonly used. The question immediately arises of how this 
extrapolation can be carried out and several methods ||21|| are employed. These are 

Pade Analysis: 

E(t) is assumed to be representable by a [L/M] Pade approximant with L = M 
i.e. it has only separable poles and zeros of finite order in any finite region of the 
complex t-plane, and the extrapolation t — > oo made directly, 

D-Pade Analysis: 

dE(t)/dt is assumed to be representable by a [L/M] Pade approximant with 
M > L + 2, on the basis of the same considerations as above, and the ground 
state energy is found by numerically integrating this 



roc 

E(oo) - E(0) = / dsE'( 
Jo 



(3) 



Laplace Method/Resolvent Analysis: 
The Resolvent 



G M/« s <T+15f>' (4) 

is assumed to have a pole structure of the form 

so that a [L/M] Pade analysis with M = L + 1 can be made and the pole found, 
Inversion Analysis: 

One inverts the E(t) function, forms dt(E)/dE and makes a Pade analysis of 
this. The ground state energy is the location of a pole in this derivative, from 

(6) 

dE E-Eq' 1 ' 

U E of F" /Partition Function: 

The function S(t) = 1 — (e~ tH ) is taken as the independent variable instead of 
t. One inverts the S(t) function and makes a diagonal Pade analysis of E{5) and 
extracts the ground state energy from 

E = lim — w . (7) 

s-*i- {■[ x\ at 

The Connected Moments Expansion, of which there are many variants, is derived 
by assuming 

n 

E{t) ~ E + • ( 8 ) 

i=i 

The variants give the following explicit forms for the ground state energy 
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CMX-HW: 

The first method gives the ground state energy up the nth order as 



Eq 



v\ - 



E 



s 2 



; 2 



1 ]lj = l Sjr 



where the S'-variables are recursively generated by 



Si,fc 



Q Q Q2 

>~>i,ki->i,k+2 ~ <Jj 



fe+1 > 



CMX-SD: 



Another methodpOfl gives the ground state energy up the nth. order as 

Eq = U)i >n , 

where the u-variables are recursively generated by 



V 3 



UJ j+ljk UJ2,k 
Uj,k — * r 



CMX-LT: 

This method described in Ref. ^] gives 



UJ i 



j+2,k 5 



V\ ~ (v 2 , V 3 , ■ ■ ■ , Vn) 



( v z 



^5 



\ ^n+l ^n+2 



^n+2 
^2n-l / 



1 / V2 \ 
\ Vn J 



(9) 



(10) 



(11) 



(12) 



(13) 



There is also another variant, called the AMX[]23"|, which is very similar to the CMX- 
HW method. The CMX, in addition to its convergence difficulties, often presents 
another, although less serious problem. Usually the cumulants are nonlinear func- 
tions of some coupling constant, and the resulting ground state energy estimate above 
truncated at some order above would be rational function of the constant. In these 
problems |24|, p5|j an unphysical pole would appear for real coupling values. This pole 
was found to be exactly cancelled |23j if one included the next order, but a next order 
is not always available. 



2 The Model 



We follow the conventions and notations of Ref. p^] as closely as possible in what 
follows. For the N site chain, with N even but not a multiple of 4, the Hamiltonian is 
taken to be 

TV 

H = E + s " s ^ > ( 14 ) 

i=l 
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where the labels . . . index the sites. Following the standard treatments a transfor- 
mation from spin operators to a mixed Fermi-Bose operator algebra is performed 

ala t = S*±tSV, (15) 

and then a Wigner- Jordan transformation to a pure Fermionic representation 



(16) 



so that the Hamiltonian is now 

H -- 



i N 
i=i 



Q+i + Cj +1 c. 



(17) 



In doing so we have neglected the boundary terms arising from the exact mapping of 
the spin or a-cyclic problem, and will only consider the c-cyclic problem Cat+i = C\. 
We take this course because we don't want to complicate matters unnecessarily and 
are only interested in the thermodynamic limits. 

Proceeding with the diagonalisation of this quadratic form via the canonical trans- 
formation 



i 

4 = E 

i 

with transformation coefficients 

4> P i ~- 



4>pi+i> P i + , 4>pi-*P P i 

n c + n — < 



(18) 



1>: 



N 

N 
cosp 



pi 



sin pi 

cos pi 

P 



p > 

P < o 



A Ypi 



(19) 



where the quasiparticle energy A p = |cosp| and the momentum values, labelled by 
p,q, . . . are 

2ixm 



p 



N 



m = -N/2, . . . , 0, 1, . . . , N/2 - 1 



The Hamiltonian then becomes 

p p 

The transformation coefficients are orthonormal in the sense 

4>p ■ 4>p' = Sppi 
i> P ■ Vy = <V 
4> P ■ Vy = o . 



(20) 



(21) 



(22) 



In all the moment formalisms a choice has to be made as what trial state the 
moments of the Hamiltonian are computed with respect to and it is this choice which 
is the crucial physical one. The only general requirement is that, for ground state 
properties at least, the trial state should have a non-zero overlap with the true ground 
state. In all the previous and approximate treatments of the the XY chain, and most 
of the other antiferromagnetic spin models, the trial state was taken to be the classical 
Neel state. This state is in fact, as we shall quantify, as far from the true ground 
state as one could possibly be but its simplicity allows for ease of moment computa- 
tion, certainly compared to alternatives. So as a "natural" choice we proceed by only 
considering this choice. 

In the spin picture the Neel state \N) will be represented by 

I UU ■ ■ • Titi) , (23) 

and in the fermionic c— picture by 

|1010 . . . 1010) . (24) 
The defining equations for the Neel state are 



Vi (25) 



4 W \N) = 
c 2i \N) = 0^ 

and using the transformation coefficients, Eq. (19), and the spectral points, Eq.( 
these can be shown to be equivalent to the following equations in terms of the quasi- 
particle operators 



( v l + Vp )\N) = 



VpEZ . (26) 



IP 'Ipj 

i 

Here the spectrum of allowed momentum values have been divided up into the set 



Z = [-7T, -tt/2) U (0, tt/2) , (27) 
and its complement Z, and there is a bijection from Z to Z defined by 

q = — q — 7r q < 
= -q + 7i q > . (28) 

If we denote the vacuum state for the quasiparticles by |0) then the Neel state can 
be expressed in terms of the quasiparticles by 

l qeZ qi<q2&Z 

■■■ + vlAvlA---vt N/2 vl /2 \\0) , (29) 
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and because of the anticommutation relations for the quasiparticle operators 

|iV)=2- JV / 4 exp (X>Jtj}J |0) 

\qeZ J 



(30) 

n^io>. 



qez 



One can easily show that this state satisfies the defining Eq. (|26|) . Basically the Neel 
state is created from the quasiparticle vacuum state by pairs of excitations whose 
momenta sum to ±7r. 

We will have need of some spectral sums defined on Z and these are 

VA^ = (2n)! N 
i 2 2n + 1 (n\) 2 

1 (3D 



§ A * sin(vr/iV) ' 

3 Generating Function 

We now show how to find the exact form for the generating function 

(et H ) = (N\e tH \N) = J2^T' ( 32 ) 

n=0 

where the moments \i n = (H n ). We will establish this result for a finite system of size 
N although ultimately the thermodynamic limit will be taken. It is straightforward to 
show a direct expression for a general moment 

fi n = 2- N / 2 lEZ + J2(Eo + 2A q ) n + ^ £ (Eo + 2A qi + 2A 92 )" 
I qez ' qi ^q 2 ez 

... + (£o + 2A ?1 + ... + 2A gjv/2 )4 . (33) 
in its symmetrical form. Then our initial form for the generating function is 



( e tff) = 2- N / 2 e Eo * J 1 + £ e 2A '* + — £ e 2t(A 91 +A 92 ) 

_ _ _ _|_ e 2t(A qi +...+Aq N/2 ) I ^ /g^N 



For simplicity let us define the following sum 



4 = ^e 2 "* A «, (35) 
qez 



and denote the general p-th term of the Eq. (34) by g p /p\ after expanding it. This 
would now be 



^ = i Yl e2tAn ■■■ e2tAqv ■ ( 36 ) 

P ' P ' qi ^...^g p ez 

Such sums over mutually distinct spectral points are recognisable as the cycle index of 
the Symmetric Group S p . Consequently we have the partition sum formula 

p\ ^ iilln^&z . . .iJp*p v ; 1 "' p ' 



(37) 



and find immediately the exponential sum, constituting the generating function, is 

\n+l 



p=0 ^' \n=l 



n 



(38) 



Each term g p /p\ is an elementary symmetric polynomial in the variables z p = e 2tAp and 
Eq. ( P?|) is just the way of re-expressing this in terms of the power sums /„ = ^2 p Zp. 

The generating function is then expressible in the following sequence of closed 
forms 

(e tH ) = 2- N ' 2 e tEo exp j i log + e 2tA «) J 
= cosh(tAg) 



nn (i+ 

k=0q£Z V 



\ t 2 A 2 

tt 2 (A; + 1/2) 2 9 



(39) 



> exp < — / dq log cosh(tA q ) 

n^oo I An J_ n 

This leads to the generating function F(t) associated with the t-expansion 

n r /2 

F(t) = — dglogcosh(£A„) , 
Jo 



(40) 



which is identical to the form of the free energy p6| for this system of noninteract- 
ing quasiparticles with zero chemical potential, with the identification |/3 — > t. The 
corresponding density fit), defined by F(t) = Nf(t) is, 



,/(*) 



n 

fc=0 



11/ t 2 

2 + 2V 1 + tt 2 (A; + 1/2) 2 



(41) 



This function is also the cumulant generating function, so from the taylor series 
expansion about t = we have 



^2n+l = 

_ 2 2n -l (2n)! 
V2n ~ An (n\) 2 



n>0 
B 2n N n > 1 



(42) 



where B 2n are the standard Bernoulli numbers. These cumulants coincide exactly 
with all the low order ones|27]] computed directly by combinatorial enumeration on a 
computer. 

There is also a simpler, more useful expression for the moments, other than that 
given by Eq. (33). Using the infinite product expansion for the generating function 
and expanding it in t 2 , one collects all factors contributing to t 2n term. By picking n 
such factors, all of which must be distinct, 



A*2ri 



1 /4\ n A* K 



V - - (43) 

(2n)\ n\\n*J (h (2k 1 + l) 2 ' ' ' (2k n +l) 2 ' 1 ) 

and using 

one can apply the cycle index formula to arrive at 

= we n£ ( H 'S~ 1) |ai|JV )'' • (45) 



4 Analyticity and Analysis 

In this part we analyse our results for the generating function and its closely related 
functions. The cumulant generating function F(t) and its derived functions have a 
finite radius of convergence about t — of \t\ < 7r/2, and which is evident from the 
infinite product form of Eq. (f£ID in having branch points at t — iri(k + 1/2) for all 
integer k. This finite radius of convergence about t — will then be a source of 
problems for extrapolation methods, even though the function itself is analytic in the 
neighbourhood of 9ft (£) > 0. 
Focussing on 

n r /2 

E(t) = / dq A 9 tanh(tA g ) , (46) 

7T Jo 

we verify the Horn-Weinstein Theorem is exactly satisfied 

N 

lim E(t) = . (47) 

Now the first excited state has an energy X\ = E + 2A q for a value of 

7T , 2 . . 
g = ± 2 (1± iV } ' (48) 

so that X% — E ~ 2tt/N and that the gap vanishes like 0(1 /N). The Horn-Weinstein 



theorem for the gap also gives zero. Finally the expression for the overlap [28 



|(iV|0)| 2 = exp - / dt[E(t)-E(oo)]} , (49) 



o 
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can also be shown to give the correct answer of 2~ N ^ 2 . The vanishing of this in the 
thermodynamic limit is another signal of the problems with using the Neel state. 

However while the specific predictions of the theorem are correct the basis for these 
are not, as will be evident in what follows. We require for this the large t development 
of the function E(t) and this can be obtained in a number of ways. Using the Mellin 
Inversion Theorems one can get the leading series only, so we adopt a direct expansion 
of the integrand. This yields 



W 2tt ^ m\T 1/2 v ;| 1 \ t, 

771=1 V ' ' 



+ 2 JL Y(- r e~ 2nt Y ( 1 /2)7n(2m + l)! ™? 1 2m _ 2+k 

n=l m=0 k=0 

Several comments are in order here. The above m sums are separately divergent sums, 
and consequently these expansions are asymptotic, but when combined the result is 
convergent. The leading terms are algebraic, and the exponential terms are subdomi- 
nant, with a numerical evaluation 

which is at complete variance with common assumptions discussed above. This is 
basically due to the vanishing spectral gaps. The specific leading exponent of 2 and not 
1 in the descending series of t is also significant and different from the Pade analysis 
assumptions. A diagonal Pade analysis of the taylor series for E(t) is not possible 
because the u 2n +x are zero, and so the t — ► oo limit either doesn't exist or is trivially 
zero. 

However a D-Pade analysis is possible as [M/M+2] approximants can be con- 
structed for E'(t). We display in Figure 1 the error for the ground state energy, 
computed using a series of D-Pade approximants for M values up to 24. It is clear 
that while there is initial convergence, this eventually ceases and no further improve- 
ment is evident. And the best result is never any better than 10~ 4 . Furthermore there 
are some gaps in our data where the Pade approximant had singularities close to the 
positive real axis which made the numerical integration algorithm unreliable. There 
are two sources of error in this analysis, the first being the error due to the finite Pade 
approximants of the integrand, and the second due to the numerical integration. We 
consider the latter error to be small (much smaller than 10~ 4 ) and relatively uniform 
across the series so that most of the variation is due to the truncation. 

The Inversion Method is also flawed here in that the leading order factor for t(E) 
has a branch cut singularity, not a pole 

/ 7T \ V2 1 

^-(s") 7W=Wr (52) 

We have done an inversion analysis, using a full range of Pade approximants [L/M] to 
Eq. (H), and present the error for the ground state energy in Figure 2 against L. There 
are some points not evident in this Figure, and that is the root of the denominator 
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polynomial closest to the exact ground state energy was not necessarily real, nor the 
one with the lowest real part. Here it is also clear that there is good initial convergence 
but the improvement ceases and in fact gets worse. There are some points much closer 
to the true answer than the rest but it is also true that there is clustering around this 
point and there is no way of uniquely distinguishing this point. Also going to higher 
order doesn't lead to any improvement. 

The Resolvent can be exactly written down as 



G{s)/s 



dx exp t — 



7T 



-x + 



tt/2 



dq log cosh (sxA q 



(53) 



in terms of a scaled inverse energy s = uN/ ir, although not explicitly evaluated. Instead 
we present a saddle point approximation to this integral for iV — > oo which yields 



G(s)/s 



N 
2^ 



tt/2 



-1/2 



dq A^sech 2 (sx A ? 









x exp < 




-x + / 






JO 



"Tt/2 



dq log cosh (s^oAq 



, (54) 



where the saddle point xo is given by the functional equation 

»7r/2 

dq A q tanh(sxoAq) . 



(55) 



The inverse energy lies in the range \s\ > 1 and as s — > — 1 then the point Xq — > oo. 
Using the asymptotic form for this integral, namely from the large t analysis Eq. (pT|), 
one finds in the neighbourhood of the ground state energy 



7T 



x 



and that the Resolvent behaves like 



G{s)/s ~ 2~ N ' 2 



4.6 3 / 2 



7T 



7/2 



A' 



2^ 



-1/2 



-1/2 



(56) 



+ 1I-3/4 



exp 



N 

2^5 



+l|l/2 



(57) 



So that we have a branch cut and essential singularity, and not simple separated poles. 

For the U E of F" analysis the inversion of the function 8(t) for 5 — > 1 _ has the 
following leading order term 



t = | log 2 + £ log(l - 6) + 0(1/ log(l - 5)) 



(58) 



and consequently the energy function E(5) has the leading term as 5 — » 1 (with fixed 
N) of 



^(5) - E 



N 3 



24tt [log(l-J)] 1 



(59) 
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which is not a rational function of 5. 

The Connected Moments Expansions actually exhibit pathological behaviour in 
this Model in that for both the CMX-HW and CMX-SD the recursion can't even get 
started because z/3 = 0. And for the CMX-LT when the dimension of the matrix to 
be inverted is odd then it is exactly singular, and when it is even the inner product is 
exactly zero, which is again due to u 3 vanishing. This behaviour is typical in that the 
expansions work better when z/ 3 is large, say relative to v^ 2 . 

5 Conclusions 

We have given the exact solution to the moment problem for the one dimensional 
isotropic spin-half antiferromagnetic XY Model and explicitly found expressions for 
the moments, the generating function and the Resolvent. With this Model we have 
shown the Theorems used by the t-Expansion and Connected Moment Expansion are 
correct but that the assumptions upon which they are based, are not true. In fact 
we find leading order algebraic decay in E(t) — E of second degree and not sums of 
exponential decay terms, that the inversion of E(t) has a leading order branch cut, and 
that the Resolvent in the neighbourhood of the ground state energy has branch cut and 
essential singularities. For this Model we find the CMX methods and some of the t- 
expansion extrapolation methods to give pathological results, and the two t-expansion 
methods that can be applied give clear indications of having asymptotic convergence 
behaviour with respect to the truncation order. The analytic properties of this Model 
are at complete variance with commonly held assumptions made by the methods and 
shows why extrapolation and expansion procedures based on these assumptions have 
poor convergence. However we optimistically predict that revision of the extrapolation 
methods, along lines indicated by our results, should dramatically improve the conver- 
gence characteristics of the t-Expansion and CMX-type methods. Although the XY 
Model is perhaps the most difficult case to treat with moment methods employing the 
Neel state as the trial state, being the most extreme quantum limit of the anisotropic 
antiferromagnetic Heisenberg Models, a number of features found here are expected 
to survive for the isotropic Heisenberg Model. In this last case the excited state gap 
still vanishes, as in our current example. We also want to emphasize that the general 
questions of convergence in these methods, which are primarily of a physical and not 
just a technical mathematical nature, still remain open. 
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2 The error for the ground state energy computed using an Inversion anal- 
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